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ON THE HARMONIC MEASURE OF STABLE PROCESSES 


CHRISTOPHE PROFETA AND THOMAS SIMON 

Abstract. Using three hypergeometric identities, we evaluate the harmonic measure of a 
finite interval and of its complementary for a strictly stable real Levy process. This gives a 
simple and unified proof of several results in the literature, old and recent. We also provide 
a full description of the corresponding Green functions. As a by-product, we compute the 
hitting probabilities of points and describe the non-negative harmonic functions for the 
stable process killed outside a finite interval. 


1. Introduction and statement of the results 

Let L = {Lt, t > 0} be a real strictly a—stable Levy process, with characteristic exponent 

^(A) = log(E[e'^^i]) = _ A e M. (1.1) 

Above, a G (0, 2] is the self-similarity parameter and p = P[Li > 0] is the positivity param¬ 
eter. Recall that when a = 2, one has p = 1/2 and d/(A) = — A^, so that L is a Brownian 
motion. When a = 1, one has p G (0,1) and L is a Cauchy process with a linear drift. When 
a G (0,1) U (1, 2) the characteristic exponent reads 

^(A) = — A|“(l — i/3 tan(7ra/2) sgn(A)), 

where (3 G [—1,1] is the asymmetry parameter, whose connection with the positivity param¬ 
eter is given by Zolotarev’s formula: 

p = -h —arctan(/3 tan(7rQ;/2)), 

2 Tia 

and = cos(7ra(/3 — 1/2)) > 0 is a scaling constant. We refer e.g. to Chapter VIII in [2] 
for more details on this parametrization. One has p G [0,1] if a < 1 and p G [1 — 1/a, l/a] if 
a > 1. When a > l,p = 1/a or a < l,p = 0, the process L has no positive jumps, whereas 
it has no negative jumps when a>l,p=l — l/aora<l,p = l. 

Set L = —L for the dual process and p = 1 — p for its positivity parameter. Throughout, it 
will be implicitly assumed that all quantities enhanced with a hat refer to the same quantities 
for the dual process, that is with p and p switched. We denote by Pj, the law of L starting 
from x G M. Introduce the harmonic measures 

Hx{.dy) = e dy, T < oo] and Hl{dy) = e dy, T* < oo]. 
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where T = inf{t > 0, \Lt\ > 1} and T* = inf{t > 0, \Lt\ < 1}. Observe that by spatial 
homogeneity and the scaling relationship 

{{kit, t >0},F^) = t>0},¥kx), k>0, (1.2) 

we can deduce from the expression of the harmonic measure of the complementary of any 
closed bounded interval, whereas the knowledge of H* gives that of the harmonic measure 
of any open bounded interval. Introduce the following notation 

= max(x,0), Can = ^ and = (t — 

TT 

In the remainder of this section it will be implicitly assumed that L has jumps of both signs. 
The corresponding results where L has one-sided jumps, which are simpler, will be briefly 
described in the last section. 

Theorem A. (a) For any x G (—1,1), the measure Hx{dy) has density 

h{x,y) = Ca,p{l + xrP{l-xrP{l+y)-^P{y-l)-^P{y-x)-^ 

ify>l and h{x,y) = h{-x, -y) if y < -1. 


(b) For any x G [—1,1]'^, the measure H*{dy) has density 
h*{x, y) = Ca,p (1 + - y)~^ - {a - 1)+^ iJa,p{t) dt 

if X > 1, and h*{x, y) = h*{—x, —y) if x < —1. 


In the symmetric case, these computations date back to [5] - see Theorems A, B and C 
therein. Notice that the results of [5], which rely on Kelvin’s transformation and the principle 
of unicity of potentials, deal with the more general rotation invariant stable processes on 
Euclidean space. In the general case. Part (a) of the above theorem was proved in Theorem 
1 of [IT], whereas Part (b) was recently obtained in Theorem 1 of [T2|. Both methods used in 
[TT] (coupled integral equations) and in [12] (Lamperti’s representation and the Wiener-Hopf 
factorization) are complicated. In this paper we show that the original method of [S] works 
in the asymmetric case as well, thanks to elementary considerations on the hypergeometric 
function 


2-^1 


a b 
c 


E 

n>0 


{a)n{b)n n 
ic)n " • 


More precisely, we use three basic identities for the latter function, due respectively to Euler, 
Pfaff and Gauss, allowing to perform a simple potential analysis of the function 


(p(t) = (1 - t)-“'’(l + t)-"/" 


(1.3) 


and to obtain the required generalization of the key Lemma 3.1 in [5]. 


Dehne next the killed potential measures 


G^{dy) = E,, 


uo 


^{Lt&dy} 


dt 


and Gl.{dy) = E^, 


r-T* 


uo 


^{Ltedy} 


dt 
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It is easy to see from the absolute continuity of the two killed semi-groups with respect 
to the original stable semigroup, that both these measures are absolutely continuous. We 
denote by g{x,y) and g*{x,y) their respective densities on (—1,1) and [—1,1]'^, the so-called 
Green functions. These functions are of central interest because they allow to invert the 
stable inhnitesimal generator on (—1,1) and on [—1, - see e.g. Formula (1.42) in [6] in the 

symmetric case. Observe that they are related to the harmonic measure and to the density 
of the Levy measure of L: 

^(y) = r(a-|-l)|l/| “ ^ {Ca,p^{y>0} + Ca,pl{y<0}) , (1-4) 

through the integral formulae 

h{x,y) = / g{x,v)iy{y — v) dv resp. h*{x,y) = / g*{x,v)v{y — v) dv 

for all X G (—1,1) and y G [—1,1]'^ resp. for all x G [—1,1]'^ and y G (—1,1), which are both 
instances of a general formula by Ikeda-Watanabe - see Theorem 1 in jS]. For this reason, 
the density of the harmonic measure coincides with that of the Poisson kernel - see M pp- 
16-17. The closed expression of the Poisson kernel and the Green function for (—1,1) in 
the symmetric case, and more generally for the open unit ball in the rotation invariant case, 
are classic results dating back to Riesz [I511I6]. We refer to [6] pp. 18-19 for more details 
and references, and to the whole monograph [6] for several extensions, all in the rotation 
invariant framework. 


Theorem B. Set z = z(x,y) = 

(a) For every x G (—1,1), one has 

g{x,y) = 


x-y 


for every x ^ y. 


P(ap)P(ap) 


y-X 


ol—1 


dt 


ifye {x, 1), and g{x,y) = g{y,x) if y e (-l,a;). 


(6) For every x > 1, one has 

2i-a 


g*ix,y) = 


T{ap)T{ap) 


{y - / i>a,pi^)dt - (a - 1)+ / 'ifa,p{t)dt / ifa^pit) dt 


ifye {x,oo),g{x,y) = g{y,x) if y ^ (l,a;), and 


g*{x,y) = 


c - 2 ^““ 


Ca,pT{ap)V{ap) 

ify< -1- 


(x-y^ / - (a - 1)+/ fJa,p{t)dt dt 


Observe that in Part (b) of the above result, the condition a: > 1 is no restriction since 
by duality we have g*{x,y) = g*{—x, —y) for every x < —1 and y G [—1,1]'^. Part (a) was 
obtained as Gorollary 4 of [5] in the symmetric case, and as Theorem 1 of ra in the general 
case. Part (b) was proved as Theorem 4 in [12], in the only cases a < 1 and x,y > 1. The 
methods of [l2l[T3|, relying on the Lamperti transformation and an analysis of the reflected 
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process, are complicated. In this paper, we observe that all formulae of Theorem B can be 
quickly obtained from the Desire Andre equation and one of the two simple lemmas leading 
to the proof of Theorem A. 


The explicit knowledge of the Green function has a number of classical consequences. In 
this paper we will focus on two of them. The hrst one deals, in the relevant case a > 1, with 
the hitting probability p{x,y) = f’xlTy < T], where Ty = inf{t > 0, = y}. 


Corollary 1. Assume a > 1 and set z = 

p{x,y) = (a - 1) 
ifx>y, and p{x, y) = p{-x, -y) if x < y 


x-y 


X — y 


. For every x,y E (—1,1), one has 

a—l nz 

/ dt 


Observe that the above formula extends by continuity on the diagonal, with the expected 
property that ^x[Tx < T] = 1. Of course, this follows from the fact that {x} is regular for x 
in the case a > 1. When a —)■ 2, Corollary 1 amounts to the very standard Brownian formula 

P,[T, <Ti] = 

1-y 

By the Markov property, one can deduce from Corollary 1 the harmonic measure of the set 
{y} U [—1,1]^. Using one of our three hypergeometric identities, it is also possible to derive 
the asymptotic behaviour of ^x[Ty < T] when x ^ y, which is fractional. Last, by spatial 
homogeneity and scaling, we can quickly recover the statement of Theorem 1.5 in [12]. See 
Remark 6 below for more detail. 


We next consider non-negative harmonic functions on (—1,1), which are the non-negative 
solutions to 

— 0 

on (—1,1), where Ca,p is the inhnitesimal generator of L. As in the Brownian case, an 
equivalent characterization - see e.g. [6] p. 20 in the symmetric case - is the mean-value 
property, which reads Kx[u{Lr,j)] = u{x) for every open set U whose closure belongs to 
(—1,1), where tu = inf{t > 0, Lt ^ U}. 

Corollary 2. The non-negative harmonic functions on (—1,1) which vanish on [—1,1]'^ are 
of the type 

X ^ X{l-x)^f’{l + x)^P-^ + /i(l + a:)“^(l 

with X, p >0. 

This result might be already known - compare e.g. with Theorem 10 p.569 in [18], al¬ 
though we could not hnd it written down explicitly in the literature. Recall that in order to 
obtain all non-negative harmonic functions on (—1,1), one needs - see e.g. Theorem 2.6 in 
[0] in the symmetric case - to add to the above functions the integral of the Poisson kernel 
h{x,y) along some suitably integrable measure on [—1,1]'^. 
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Both Corollaries 1 and 2 could be obtained for the process killed inside the interval (—1,1), 
with analogous computations relying on Part (b) of Theorem B. But the formulae have a 
rather lengthy aspect, so that we prefer leaving them to the interested reader. The remainder 
of the paper is as follows. In the three next sections we prove Theorem A, Theorem B, and the 
two Corollaries. In the last section we gather, for the sake of completeness, the corresponding 
formulae in the cases of semi-hnite intervals and of one-sided jumps. 

2. Proof of Theorem A 

As mentioned in the introduction, the argument hinges upon three classical hypergeometric 
identities, to be found in Theorem 2.2.1, Formula (2.2.6) and Formula (2.3.12) of pQ, and 
which will be henceforth referred to as Euler, Pfaff and Gaus£| formula respectively. 

2.1. Proof of Part (b). 

2.1.1. The case a < 1. We reason along the same lines as in Theorem B of [5]. Set Pt{x) for 
the transition density of L. The following computation, left to the reader, is a well-known 
consequence of fll.ll) . Fourier inversion and the Fresnel integral: one has 

poo 

/ pt{z)dt = - a)ca,pZ°‘~^ 

Jo 

for every ^ > 0. Observe that by duality, one also has 

poo poo 

/ Pt{z)dt = / Pt{-z)dt = T{1 - a)ca,p\z\'^~^ 

Jo Jo 

for every z < 0. Applying the Desire Andre equation (2.1) in [5] and letting s ^ 0 therein 
shows that 


u{t,y)H*{dt) = Ca,p\x-y\ 


a—1 


( 2 . 1 ) 


'-1 


la—1 


for every x > 1 and i/ G (—1,1), where we have set 

In the symmetric case, this Abelian integral equation with constant boundary terms is solved 
in Section 3 of |5], following the method of [15]. See also [^ for the original solution, with a 
more general term on the left-hand side. After proving the following lemma, which remains 
valid for a G (1, 2), we will see that the pole-seeking method of [T5] applies in the asymmetric 
case as well. 

Lemma 1. The unique positive measure on (—1,1) satisfying 


u(t,y)ix(dt) = 1, 9 e (-1,1) 


( 2 . 2 ) 


'-1 


has the density pit) given in fll.3p . 


^among of course many others. This one is a simple consequence of the two-dimensional structure of the 
space of solutions to the hypergeometric equation. Notice that it can also be obtained by Mellin-Barnes 
inversion. See the end of the article Calculs asymptotiques in Encyclopedia Universalis. 
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Proof. The fact that there is a unique measure solution of fl2.2p is a standard fact in potential 
theory - see e.g. Theorem 1 in [13] or Proposition VI. 1.15 in [3]. In our concrete context, 
this unicity can also be obtained by a straightforward adaptation of Lemma 4.1 in [5]. To 
show the lemma, we compute by a change of variable 

/•l ^ r^+y 

/ + = Ca,p 


'-1 


-l-y 


+ C, 


a,p 


f""Xl -y- + y + t)-^P dt. 


Using two further changes of variable, the Euler formula, and the complement formula for 
the Gamma function, we transform the expression on the right-hand side into 


r(a) 


1 + 9 


ap 


r(ap)r(l-|-ap) \1 — y^ 

and then, using the notation 


iFi 


ap a y + 1 
1 -|- ap’ y — 1 

y + 1 

^ = -h’ 

y- 1 


+ 


p f^-y 


p\^ + y 




ap a p — 1 
_1 -1- ap y + l 


into 


T{a){—z)°'P 



a ap 
, ; ^ 

1 P(-^) “ p 

+ . 2-^1 

a ap 1 

1 -|- ap 

p 

1 -|- ap 2 ; 


r(ap)r(l -1- ap) 

where the last equality follows from the Gauss formula. 


= 1, 


□ 


Remark 1. The solution to (12. 2 h in the symmetric case was obtained in Lemma 3.1 of [5], 
via a reflection argument. Alternatively, the non-symmetric solution can be deduced in a 
constructive way, following the approach of [I5] pp. 41-42 or that of [7]. Observe that the 
above argument is signihcantly shorter than in these three references. 


We can now hnish the proof. Introduce the changes of variables 


t = X + 


X 


and 


y = X + 


X 


(2.3) 


X — S X — z 

and observe that they map (—1,1) onto (—1,1), in a decreasing way. Plugging these changes 
of variables into fl2.2l) implies after some computation that 

A-i 


I a—1 


(x + l)“^(x - 1)“^ J {1 +s)-'^P{x-s)-^u{s,z)ds = jx-z 

for every x > 1 and 2 ; G (—1,1). Multiplying both sides by Ca,p shows the required solution 
to fl2.1l) . which is unique by Lemma 1 and the changes of variables fl2.3p . 

□ 

Remark 2. In the following, we shall make a repeated use of the changes of variable fl2.3p . 
which may be written formally : 


|1 + X 


ap I 


X 


ap 


\y-t\ 


,_iii11"“^ 

\x —1\ 


dt = \x — y\ 


a—1 


I a—1 


(p{s)ds. 
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The interest of this change of variable is to transform an Abelian integral with two inside 
parameters into an integral of the hypergeometric type, with one parameter inside. 


2.1.2. The case a > 1. We follow the method of Theorem C in [5]. Recall that since L 
a.s. hits points in hnite time, the measure H*{dt) has total mass one. We will need the 
evaluation 

poo 

/ (Ptiz) - PtiO)) dt = r(l - a) 


for every z > 0, which is easy and classical - see the introduction of 


This implies 


e ^pt{z)dt — pi(0)r(l — J T{l — a)ca,pZ 


as s ^ 0, for every z > 0. Proceeding as in [5] pp. 544-545 shows that 


1 


Ca,p\x-y\'^ ^ = 
for every x > 1 and y G (—1,1), where 


'-1 




(2.4) 






is a non-negative function which will be determined in the same way as in (4.1) of [5]. 
Multiplying both sides of fl2.4p by (p{y) and integrating on (—1,1) shows by Lemma 1 that 

I T{y)dy'^ {ca,p j {x-y)'^-^p{y)dy - 1 

for every x > 1. The next lemma, generalizing the second part of Lemma 3.1 in |5j, allows 
to compute the right-hand side. 




Lemma 2. With the above notation, one has 




{x-yY ^{y)dy = 1 


r(l - ap) 2 ^-" 


ipa,p{t) dt 


J-I r(ap)r(l-a) 

for every x > 1. 

Proof. As in Lemma 1, a change of variable and the Euler formula show hrst that 


• / f \ ^a-l / 4 .7 r(l-ap) /x + 1 

sm Trap / {x - y) p{y) dy = -: 

T{ap)T{2-a 


a—1 


2F1 


r(l — ap) 


X 


a—1 


2F1 


1 — a 1 — ap 2 
2 — a ’ X -|- 1 
1 — a 1 — ap 2 
2 — a ’ 1 — X 


r(ap)r(2 — a) \ 2 

where the second equality follows from the Pfaff formula. Using now the Gauss formula, we 
next transform 


[-zY-YF, 


1 — a 1 — ap 1 
2 — a ’ z 


P(ap)P(2-a) ^ (ffJzJlf-zYYFi 


P(1 — ap) 


ap 


— ap ap 
1 + ap 


-,z 
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with the notation z = {1 — x)/2. Putting everything together and applying again the Euler 
formula completes the proof. 

□ 

We can now hnish the proof of the case a > 1. From Lemma 2 and an easy computation, 
we hrst deduce 


1 


1 — ap ap 1 — X 
1 + ap ’ 2 


Coming back to (12.dh and reasoning as in [5] p.552, we finally see from Lemma 1 that H*[dy) 
has density 

To conclude the proof, it suffices to observe by the Euler formula and a change of variable 
that 


1 f X — 1 


= c„,p(l + |/) “'’(1-1/) “'’(a-1)2“ ^ 


= Ca,p{l + y) “^(l-i/) “^(a-1) / 'ijja,p{t) dt. 


ap 


2 F 1 


1 — ap ap 1 — X 

1+ ap ’ 2 


□ 


Remark 3. Since k* ^(a;) is finite and positive, we can deduce from Karamata’s Tauberian 
theorem that 




r(l — a) sin(7r/a) 
7rpi(0) 




)tV 


Q—1 


as t —)■ +CX 0 . 


This asymptotic is given in Corollary 3 of [S] in the symmetric case, and in Theorem 2 of 
[IT] in the asymmetric case, with a more general formulation. Notice that T* has inhnite 
expectation. 


2.1.3. The case a = 1. This case is known to be more subtle from the computational point 
of view, because it involves logarithmic kernels. The transition density of Lt is 

= ,2 I 9. ^ 

F + 2tx cos up + x^ 

The process L does not hit points a.s. but it is recurrent, so that H*{dt) has total mass one. 
After some computation, one hnds 

- pt{x)) dt = Ci,plog|a:|. 

See also na p.391. With this formula, it is possible to finish the proof as in the case a > 1, 
but the computations are lenghty and we hence prefer to invoke a simple argument relying 
on the Skorokhod topology. Fix p G (—1,1) and let a 1. It follows from Corollary VIL3.6 
in [To] that 



C{L^'P) => C{L^'P) 
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with obvious notation for and and where ^ means weak convergence in the classical 
Skorokhod space. Using Remark VI.3.8 and Proposition VI.2.12 in (TU], it is then easy to 
deduce that 


r“>P 

±Jrp 



The conclusion follows from pointwise convergence of the densities h*{x,y) as a 1, and 
Scheffe’s lemma. 

□ 


Remark 4. The above argument relying on a.s. continuity for the Skorokhod topology will 
be used repeatedly in the sequel, under the denomination Skorokhod continuity argument. 


2.2. Proof of Part (a). By the Skorokhod continuity argument, it is enough to consider 
the case a 1. Fixing x G (—1,1) and proceeding as in [5] pp. 544-545, the harmonic 
measure Hx{dt) is seen to be the unique solution of the equation 


u{x,y) 


'(-1,1)'= 


u{t,y)H^{dt) 


(2.5) 


for every y G [—1,1]'^. In the case a < 1, this is indeed an immediate consequence of the 
Markov property, leading to the corresponding equation fl2.ip . And in the case a > 1, the 
well-known fact - see Lemma 4.1 in [I9] - that the tail distribution of T is exponentially small 
at inhnity implies that the perturbative term Ka,p is zero in the corresponding equation fl2.4p . 
Dehne 


fix {,dt^ 


Ca,p^{t)dt iff < X, 

Ca,p 0ii) dt if t > X. 


We shall deal with the two cases y > 1 and y < —1 separately. 


(i) Let 1 / G (—l,a;). Applying Lemma 1 with p and p interchanged, we get 


^1/ 


•'OL.p 


H^x{.dt) + Ca,p / y^xidt) + Ca,p / ^xidt) 


'-1 
= C. 


a,p 


c«,pk-^r 0{t)dt + / c«,p|z/-t|“ ^(p{t)dt] = Ca,p. 




The changes of variable fl2.3p implies after some rearrangement 



^Hx{dt) 


Ca,p{y-xY ^ 


for every y > 1 with the required expression for Hx{dt), which is Equation fl2.5p . 
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(ii) Take now u G {x, 1). Applying again Lemma 1, we have 


nx 


-'a^p 


y - + c, 


a,p 


y - + c, 


a,p 


y-t\°' iix{dt) 


'-1 
= c, 




-o,p| 


t|“ !p(t)dt + / Ca,p\y —1\'^ (p{t)dt 




= C, 


OL^P' 


The same changes of variables (12.3p gives 


'[- 1 , 1 ]= 


Hx{dt) = Ca,p{x-y) 


OL—1 


for every i/ < — 1, which is again Eqnation (12.51) . 


□ 


Remark 5. (a) The behavionr at inhnity of the distribution function of T is more mysterious 
than that of T*. In the non-subordinator case it is known - see Proposition VIII.3 in [2] - 
that there exists positive and finite such that 

-logP^[T > t] ~ —Kxt as t ^ + 00 , 

but the exact value of k,x is unknown except in the completely asymmetric case - see [3]. We 
refer to Chapter 4 in [6] for more on this topic in the rotation invariant case. Notice that 
the result of Theorem B (a) allows to compute the expectation of T: 




r(Q' + 1) 


(b) With our computations, we can also check the values of the total masses Hx{—1, 1)^ 
and Tr*(—1,1). On the one hand, Lemma 1 and the change of variables fl2.3l) imply 

Hx{dt) = Ca^p / {x — z)‘^~^ip{z)dz + Ca,p / {z — x)°‘~^ip{z)dz = 1. 


'(- 1 , 1 )= 


'-1 


On the other hand, in the case a > 1, fl2.3p and Lemma 2 show that 


'-1 


H*{dt) = Cc,,p / ip{z)dz - <p(x) / ip{y)dy 


'-1 


'-1 


= 1 


because 


( r(lp)raL) + (»)%)/ = i. 

2i-«r(l - ap) 


(p{y)dy = 2^ " B(1 - ap, 1 - ap) = 


j-i (1 - a)c«,pr(ap)r(l - a) 

In the case a = 1, the measure H* has also total mass one by continuity. In the case a < 1, 
we hnd 

r H*{dt) = 1 - P,[T* = cx)] = 1 " 


'-1 


r(ap)r(l — a) 


pa,p{t)dt, 
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in accordance with Corollary 2 of [5] and Corollary 1.2 of |12] . 

3. Proof of Theorem B 

3.1. Proof of Part (a). It is enough to consider the case y > x, the case x > y following 
from Hunt’s switching identity - see e.g. Theorem II.5 in [2]. By the Skorokhod continuity 
argument, it is also enough to consider a ^ 1. Reasoning as above, the Desire Andre equation 
yields 



(^Ca,p{y - x) 


L 


-1 



OL—l 


—oo 


with Ca = r(l — a). Changing the variable as in (12.3p . we deduce 


9{.x,y) 



and the result follows from Lemma 2, since z > 1. 


□ 


3.2. Proof of Part (b) in the case a < 1. 


3.2.1. The case y > 1. Hunt’s switching identity shows again that it is enough to consider 


the case y > x. As above, the Desire Andre equation and the change of variable fl2.3p imply 


9*ix,y) 



with z > X > 1, and we can conclude by Lemma [2l 


3.2.2. The case y < —1. Still using fl2.3p . we now have 


9*{.x,y) 



r(l - a) (x -yY ^ (^1 - j 





with e (1,3;), and we again conclude by Lemma [2l 


□ 
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3.3. Proof of Part (b) in the case a > 1. We only consider the case y > x. The case 
X > y > 1 is obtained by Hunt’s switching identity and the case y < —1 by analogous 
computations. Proceeding as for Equation (12.dh . we hrst deduce 

g*{x,y) = r{l-a)co,,p (^{y - x^-^ - j {y - - r(l - a)<^^(x). 

Using Theorem 1 and the computations of the case a < 1, the expression may be transformed 
into 

r(^p)r(«p) (^) / - pi (i - . 

The result follows from the hat version of Lemma 2 and the expression of K'^pix). 

□ 


4. Proof of the Corollaries 

4.1. Proof of Corollary 1. By duality, it is enough to consider the case x > y. From Part 
(a) of Theorem B and a change of variable, we see that g{x, y) extends by continuity on the 
diagonal, with 

^ (a - l)P(ap)P(ap) 

Moreover, it is clear that g vanishes on the boundary {|a:| = 1} U {|i/| = 1} and is hence 
bounded on (—1,1) x (—1,1). By Proposition VI.4.11, Exercise VI.4.18 and Formula V.3.16 
in [1], we deduce 

9{x,y) 

9{y.y) 

and the conclusion follows by Theorem B. 

□ 

Remark 6. (a) In the case « < 1, the process L does not hit points, so that the problem 
is irrelevant. In general, one can ask for an evaluation of the probability < T] where 

/ is a closed subinterval of (—1,1) not containing x, and Tj is its hrst hitting time. In the 
transient case a < 1, this problem is solved theoretically as a particular instance of the 
so-called condenser problem - see Formula (3.4) in [8]. It is an interesting open problem to 
hnd out an explicit formula in the real stable framework. 

(b) By the Markov property, we can write down the following expression for the harmonic 
measure Hx^\dt) of the set {y} U [—1,1]'’ : 

H^y\dt) = p{x,y){6iy}{dt) - Hy{dt)) + Hx{dt). 

In particular, for every x,y G (—1,1), one has 

Fx[Lt e dt,T < Ty] = Hx{dt) - p{x,y)Hy{dt). 


i?x[Ty<T] 



(4.1) 
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(c) It is interesting to mention that using the Gauss formula, we can deduce the asymptotic 
behaviour of > T] when x —)■ y, which is fractional. For instance, if i/ = 0, one has 


P,[To>T] - 

x—>-0+ 


r(2 — a)r{ap) 
r(l — ap) 


(2x 


iQ —1 


and ¥^[To>T] ~ 

x^-O- 


r (2 — a)T{ap) 
r(l — ap) 


|2x| 


OL—1 


(d) By fll.2l) and spatial homogeneity, it is easy to deduce from Corollary 1 the following 
expression of p(x, y) = ¥x\Ty < r] where r = inf{f > 0, > 1} : one hnds 


p{.x,y) 


{a 


1 ) 


x-y 

1 - 2 / 



1 — x I 


a X > y, and p{x,y) = p{—x,—y) ii x < y. When y = 0, this is Theorem 1.5 in [T2] . 
correcting a misprint (the 1 — 1/x in the second integral should be — 1/x) therein. Notice 
that Corollary 1.6 in [12] is also analogously recovered from fl4.ip . 


4.2. Proof of Corollary 2. By the general theory of Martin boundary - see e.g. Theorem 
1 in m, we need to compute the Martin kernels 


Ml (x) = lim 

y^l 


9{x,y) 


and M_i(x) = lim 


9{x,y) 


9{0,y) 


Part (a) of Theorem B and a straightforward asymptotic analysis show that these Martin 
kernels exist and equal respectively 

Mi(x) = (1-x)“^-^(l+ x)“^ and M_i(x) = (1 + x)“^-^(l - x)“^ 


whence the result. 


□ 


5. Final remarks 


In this section, we briefly describe the analogues of the above results in the case of semi- 
hnite intervals and in the spectrally one-sided situation. 


5.1. The case of semi-finite intervals. By scaling and spatial homogeneity, one can 
deduce from Theorem A - either its Part (a) or its Part (b) - the following expression of 
the density of Lr under P^,, where x < 1 and r = infjf > 0, Li > 1}. One hnds 


hriy) 


Ca,p{l - x)°‘P 

{y - 1 )"^( 2 / - x) 


This expression has been found by several authors and can be obtained in different ways 
(see Exercise VIII.3 in [2] and the references therein). Observe that it serves as a starting 
formula in na in order to prove Part (a) of Theorem A. Notice last that the expression 
extends to the case with no negative jumps, by the Skorokhod continuity argument. In the 
relevant case with no positive jumps a > 1, p = l/a, the law of is a Dirac mass at one. 

The Green function is 


9T{x,y) 


(y-x) 


a—1 


P(«p)P(a/5) 


1 -?/ 
y — x 


'lpa,p(t) dt 
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ii X < y <1 and Qrix, y) = griy, x) if 1 / < a; < 1. In the case a > 1, the analogue of Corollary 
1 which is already given in Remark 6 (d) above, can then be recovered. Finally, one hnds 
that the non-negative harmonic functions vanishing on (1, -|-oo) are of the type 

A(1 

with A,/i > 0, in accordance with Theorem 4 in [18] and the paragraph thereafter. 

5.2. The case of stable processes with one-sided jumps. By duality, it is enough to 
consider the two cases a < 1, p = 1 and a > 1, p = 1/a. 


5.2.1. The case a < 1, p = 1. It follows readily from the above paragraph that 

for all X G (—1,1). See also Example 3 in [9] and the references therein for the expression of 
the density of {Lt-,Lt) under P^,. Similarly, one has 

Ca,l\l + X\ 


h*{x,y) = 


{l + y)^iy-x) 


for all X < —1 and h*{x,y) = 0 for all x > 1. In accordance with the fact that L is a 
subordinator, the Green function is 


( 2 /-a;)" ^ , 
9\x,y) — r(Q!) 


for all X, 1 / G (—1,1), 


9*{x,y) = 


{y-x) 


OL—1 


r(a) 




|l-|-a;|(y-l) 
2 


(1 -I- f) dt ] 1 


R>i} 


for all X < —1, and g*{x,y) = g{x,y) for all x > 1. The problem of Corollary 1 is irrelevant. 
Finally, the non-negative harmonic functions on (—1,1) vanishing on [—1,1]^ are constant 
multiples of (1 — x)““^. 

5.2.2. The case a > 1, p = 1/a. Using Skorokhod continuity in Theorem A (a) and the 
absence of positive jumps, one has 

H.{dy) = c„,i_i/,(l-x)(l+x)“-i(l- 2 /)-i|i/+l|i-“(x-i/)-'l(,<_i}di/ + P,[Ti <T]<5i(dp). 

Either taking the limit in Remark 6 (d) or integrating the hrst term, we can compute the 
weight of the Dirac mass, and hnd 

( I 1 \ a—1 

j Si{dy). 

The corresponding Green function is 

a—1 


9ix,y) = 


F(a) 


(1 -|/)(l + a;) 


{x-yr-H 


{x>y} 
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The hitting probabilities are 


P^IL<r| = (^) 


Q—1 


for every x < y, which is also a consequence of a well-known result on scale functions - see 
e.g. Theorem VII.8 in [2], and 

a —1 / c\/ \ \ 0.-1 


= (Ff) - {¥t)' 


for every x > y. Finally, the non-negative harmonic functions on (—1,1) which vanish on 
[—1, are of the type A(1 — x)^~^{l + -|- /i(l -|- with A, p > 0. 


It is clear that H*{dy) = 6_i{dy) for all x < —1. To compute H*{dy) for x > 1, let 
us introduce r* = inf{t > 0, < 1}. The absence of positive jumps and the formula for 

semi-hnite intervals imply after some computation 

H*{dy) = e dy] + < -l]6_i{dy) 

( x — 1 

- zY~‘^dz 

in accordance with Remark 3 in [T3] - see also Proposition 1.3 in [T2] . 
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\S-i{dy) , 


^a,l—l/c 


(x- 1)"-^! 


x-y 
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